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I. INTRODUCTION 



The modification of the QCD vacuum, and of its thermal excitations as well, under the influence of external fields, 
, is an attractive topic. Firstly, it is extremely interesting to understand how an external field can modify the main 
• characteristics of confinement and spontaneous chiral symmetry breaking. Secondly, strong magnetic fields, with order 
of magnitude between eB ps to J and eB 15to^, might be produced in the very first moments of the non-central 
heavy ion collisions [HH)]. In this case, it has been argued that the non-trivial topological structure of thermal QCD 
gives rise to Chiral Magnetic Effect (CME) [U, 13] ■ 

An useful approach to the physics of strong interactions in external magnetic fields is the use of some model. 
Among them, the Nambu-Jona Lasinio (NJL) model @ is quite popular, see Refs. @ for reviews. In this model, 
the QCD gluon-mediated interactions are replaced by effective interactions among quarks, which are built in order 
to respect the global symmetries of QCD. On the other hand, the NJL model lacks confinement of color. It is well 
known that color confinement can be described in terms of the center symmetry of the color gauge group and of the 
Polyakov loop Q , which is an order parameter for the center symmetry. Motivated by this property, the Polyakov loop 
extended Nambu-Jona Lasinio model (P-N JL model) has been introduced @, , in which the concept of statistical 
confinement replaces that of the true confinement of QCD, and an effective interaction among the chiral condensate 
and the Polyakov loop is achieved by a covariant coupling of quarks with a background t emp oral gluon field. In the 
literature, there are several studies about various aspects of the P-N JL model [HMLfjl [lJH26[ . Lattice studies on the 
response of the QCD ground state to external magnetic and chromomagnetic fields can be found in [27l - l3ll ] . Previous 
studies of QCD in magnetic fields, and of QCD-like theories as well, can be found in Refs. [33 - l38j . 

Beside the Polyakov loop, it has been suggested ,39j that another observable which is an order parameter for the 
center symmetry, hence for confinement, is the dressed Polyakov loop. The dressed Polyakov loop has been computed 
in Refs. (4fj| | within the scheme of truncated Schwinger-Dyson equations; within the Nambu-Jona Lasinio model, Si 
has been computed at finite temperature and chemical potential in [4lJ. Finally, the dressed Polyakov loop has been 
computed within the PNJL model in [2l| at finite temperature, and in |26l| at finite temperature with strong magnetic 
field. 

In this talk, I present results obtained in Ref. [26| about the phase structure and the dressed Polyakov loop of hot 
two massive flavor quark matter at zero chemical potential, in an external magnetic field. To compute the effective 
potential, I rely on the PNJL model of strongly interacting quarks. 

II. THEORETICAL FRAMEWORK 

I consider here two flavor quark matter whose Lagrangian density is specified as [26[ 

C = q^Df, - too) q + 9a [{qq) 2 + {qil^rq) 2 ] 

+gs[(qq) 2 + (qij5Tq) 2 ] 2 . (1) 

The covariant derivative embeds the quark coupling to the external magnetic field and to the background gluon field 
as well, see below. In Eq. (|TJ) , q represents a quark field in the fundamental representation of color and flavor (indices 
are suppressed for notational simplicity); r is a vector of Pauli matrices in flavor space; toq is the bare quark mass, 
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which is fixed to reproduce the pion mass in the vacuum, m x = 139 McV. The model at hand is called Polyakov 
loop extended Nambu- Jona Lasinio model (PN JL in the following) , since a coupling of the chiral condensate and the 
Polyakov loop is introduced via the covariant derivative in Eq. ([T]), see below. 

In this study, I limit myself to the one-loop approximation for the partition function. In order to couple the 
Polyakov loop to the quark fields, it is customary, in the PNJL model, to introduce a background temporal, static and 
homogeneous Euclidean gluon field, A4, in terms of which the Polyakov loop is given by P = Tr[exp(i/3A4)]/3. A4 is 
coupled to the quarks via the covariant derivative, see Eq. (1); as a consequence, a coupling among the quark fields 
and the Polyakov loop arises naturally when the integration over fermion fields in the partition function is performed. 

I work in the Landau gauge, and take the magnetic field homogeneous, static and aligned with the positive z— axe. 
Moreover, I take twisted fermion boundary conditions along the compact temporal direction, 

q(x,P) =e- iv q(x,0) , tp € [0, 2ir] , (2) 

while for spatial directions the usual periodic boundary condition is taken. The one-loop thermodynamic potential in 
the general case of twisted boundary conditions is given by [26| 



Q=U(P,P,T) + — + ——- ^ —^^ ak —gA(j> z ,k)u k (p z ) 
9a 9 ° f=u,d Z7r k J -°° Z7f 

-T J2 E ak / + °° — l0 S i 1 + ^Pe-^- + 3Per^ £ - + e~^ £ -) 

_ T ]T E f +Q ° l0 S C 1 + 3Pe ^ £+ + 3Pe- 2p£ + + e- 3 ^) . (3) 

f=u,d k J -°° 

In the previous equation, a = g a (qq) = 2g,j(uu); k is a non- negative integer which labels the Landau level; 
«fc = <5fco + 2(1 — Sko) counts the degeneracy of the k— th Landau level. Moreover, 0Jk(p z ) 2 = p 2 . + 2\qjeB\k + M 2 , with 
M = mo — 2a — 4cr 3 <?8 / g 3 ■ The arguments of the thermal exponentials are defined as £± = ui^ (p z ) ± i(ip — ir) / '/?, with 
f defined in Eq. The vacuum part of the thermodynamic potential, f2(T = 0), is ultraviolet divergent. In this 
study, I use a smooth regularization procedure by introducing a form factor <?aQ?) in the diverging zero-point energy. 
The potential term U[P,P,T] in Eq. ^ is built by hand in order to reproduce the pure gluonic lattice data [lOl ] . 
Among several different potential choices I adopt the logarithmic form of [Io| . 

Following [39j], I introduce the dual quark condensate, 

E n (m,V) = J o £^(qq) G , (4) 

where n is an integer. The expectation value (^q denotes the path integral over gauge field configurations. The case 
n = 1 is called the dressed Polyakov loop. For my later convenience, I scale the definition of the dressed Polyakov 
loop in Eq. (j4]), and introduce 

^ = -2ng a J 2 J ^e-^(qq) G , 

\ e-*^) , (5) 







where cr(ip) corresponds to the expectation value of the a field computed keeping twisted boundary conditions for 
fermions. 



III. RESULTS 



In the left panel of Fig. [I] I collect the results for the dressed Polyakov loop as a function of temperature and 
magnetic field strength. In the right panel of the same figure, I plot the data of the effective susceptibility, dUi/dT. 
The bifurcation of the effective susceptibility at large value of the magnetic field strength arises because the dressed 
Polyakov loop is capable to feel both the deconfinement and the chiral crossovers, which are split of almost 15% if the 
magnetic field strength is large enough. 

In Figure [2j I collect the results on the pseudo-critical temperatures for chiral and Polyakov loop crossovers, in the 
form of phase diagrams in the eB — T plane. The dashed line denotes the Polyakov loop crossover, and the dot-dashed 
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FIG. 1. Dressed Polyakov loop (left panel) and its effective susceptibility (right panel) as a function of temperature and 
magnetic field strength. 
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FIG. 2. Phase diagram of the PNJL model in magnetic field. Dashed line denotes the Polyakov loop crossover; dot-dashed line 
corresponds to the chiral crossover. The shaded area is the region, in the eB — T plane, in which quark matter is not statistically 
confined, but chiral symmetry is still broken by the chiral condensate. Temperatures on the vertical axes are measured in units 
of the pseudo-critical temperature at zero field, which is To = 175 MeV. 
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line corresponds to the chiral crossover. The shaded area is the region, in the eB — T plane, in which quark matter is 
not statistically confined, but chiral symmetry is still broken by the chiral condensate. Temperature on the vertical 
axes are measured in units of the pseudo-critical temperature at zero field, which is To = 175 MeV. Wc fit our data 
on the pseudo-critical temperatures by the law 

T C A feB\ a 

where A = a, P, and To = 175 MeV. The numerical values for T x are a = 2.4 x 10~ 3 , a = 1.85; for Tp the numerical 
coefficients arc a = 2.1 x 10~ 3 , a = 1.41. 

It is instructive to compare this result with those obtained in a different model. The shape of the phase diagram 
drawn in Fig. [5] is similar to that drawn by the Polyakov extended quark- meson model, see e.g. Fig. 13 of Ref. j37j . 
In that reference, an interpretation of the split in terms of the interplay among vacuum and thermal contribution, is 
given. I totally agree with those arguments, which are reproduced within the PNJL model as well, as the results on 
critical temperatures show. 



IV. CONCLUSIONS 



In this talk, I have reported about the computation of the dressed Polyakov loop and of the phase diagram of hot 
quark matter under the influence of a strong external magnetic field. The results on the dressed Polyakov loop, £i, 
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in magnetic field show that this quantity is capable to describe both Polyakov loop and chiral crossovers. This is 
resumed in the double peak structure of its effective susceptibility. We measure an increase of both deconfmcmcnt 
and chiral crossovers; the tiny split of the two critical temperatures is of the order of 10% for the largest value of the 
magnetic field strength considered here. 

I acknowledge the organizers of QCD@WorklO for the opportunity to give a talk. I acknowledge R. Gatto for his 
fruitful collaboration in Ref. [2^| . on which this talk is mainly based. Moreover, I acknowledge correspondence with 
M. Huang and S. Nicotri, and in particular with M. d'Elia. Furthermore, I acknowledge stimulating discussions with 
L. Campanelli and K. Fukushima. This work is supported by JSPS under the contract number P09028. The numerical 
calculations were carried out on Altix3700 BX2 at YITP in Kyoto University. 
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